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Q I Abstract 

H ' 

Context. X-ray surveys allow to identify young, main-sequence stars in the solar neighborhood. Young, stellar samples, selected according to 
their activity, can be used to determine the stellar birthrate in the last billion years. The ROSAT North Ecliptic Pole survey (NEP), with its 

1 1 1 moderately deep sensitivity (fluxes « lO"''' erg cm"- sec"'), is the best survey, to date, able to sample the intermediate-age (10** - 10' years) 
' nearby population. The identification process of NEP X-ray sources resulted in 144 X-ray sources having a normal stellar counterpart, with an 

excess of yellow stars with respect to model predictions. 

Aims. We want to determine if these X-ray active stars are young or intermediate-age stars, or active binaries. 

Methods. We acquired high-resolution, high signal-to-noise ratio optical spectra, to determine the age and physical properties of the NEP 
X-ray-detected stellar sources. We measure the (i) lithium abundance using the Lii 6707.8 A line, which is an excellent, youth indicator for 
our age range of interest; (ii) rotational and radial velocities (through cross-correlation methods); and (iii) chromospheric emission (from Ho- 
and Nai Di and D2 lines). 

Results. The radial velocities distribution is consistent with that of a young field star population of age 4x10' yrs, or younger Rotational velocity 
measurements imply that our sample is dominated by relatively young or intermediate-age stars, as confirmed by our lithium measurements. 
^ Conclusions. Most of the detected stars probably belong to a young or intermediate-age population. Our measurements suggest that a burst in 
the stellar birthrate of a factor of four occurred in the last 10** years. We cannot, however, exclude the possibility that a small fraction of sources, 
amongst the fastest of the K-rotators, are old binary systems with tidally-locked rotation. 
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1. Introduction 

The analysis of the stellar component of flux-limited X-ray sur- 
veys is a powerful tool in the study of the young stellar popu- 
lation of the Galaxy. Optical data, which helps determine pho- 
tospheric parameters, cannot alone discriminate between main- 
sequence stars of different age. In contrast, the average X-ray 
luminosity of normal, late-type stars decreases by up to 3 or- 
ders of magnitude mainly during the main-sequence lifetime 
(Micela et al. 1985, 1988; Guedel et al. 1997; Feigelson & 
Montmerle 1999); optical properties, however, remain almost 
unchanged. Stellar X-ray surveys aUow to study in detail the 
global properties of stellar populations in the solar neighbor- 
hood, in particular of the young population, since young stars 
are detected, in X-rays, out to larger distances than older stars. 
Shallow stellar, flux-limited. X-ray selected samples will be 
dominated by young stars; in contrast, deep high-latitude stellar 
X-ray samples will be dominated by old stars, because young 
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stars have smaller scale height, and deep X-ray surveys probe 
distances beyond the scale height of the youngest population 
(Micela et al. 2001). The comparison between observations 
and predictions from galactic models offers an important tool 
to constrain the spatial distribution of young stellar population 
and to infer the local stellar birthrate in the last billion years. 

Analyses of both the Einstein Extended Medium Sensitivity 
Survey (EMSS, Gioia et al. 1990), and XMM-Newton Bright 
Serendipitous Survey (XBSS, Delia Ceca et al. 2004) measured 
an excess of yellow stars (Sciortino et al. 1995; Lopez-Santiago 
et al. 2006, respectively). For the EMSS, these yellow stars 
appear to be young, in agreement with their measurement of 
lithium abundance (Favata et al. 1993). For the XBSS, the ex- 
cess may represent a young population, or a number of active 
binary systems that have a yellow dwarf as a primary (Lopez- 
Santiago et al. 2007). 

The study of the stellar content of the ROSAT North 
Ecliptic Pole survey (NEP hereafter, Henry et al. 2001) con- 
firmed the existence of an excess of yellow stars in the solar 
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neighborhood even at intermediate X-ray flux, which is proba- 
bly due to a young steUar population (Micela et al. 2007, here- 
after referred to as Paper i). 

The NEP survey is a survey of X-ray sources, covering 
a 9° X 9° area centred on the Galactic North Ecliptic Pole 
(I = 96°.4, b = 29°.8). This is the sky region with the high- 
est sensitivity observed during the ROSAT All Sky Survey 
(RASS), with a total exposure time amounting to ~ 40 ks at 
the pole. The relatively high galactic latitude, together with the 
moderate sensitivity, permits to observe young stars close to 
their scale height. 

An overview of the NEP survey, including both its selec- 
tion function and optical identification program, can be found 
in Henry et al. (2001); Voges et al. (2001) summarized the X- 
ray data and source statistical properties; Gioia et al. (2001) 
provided evidence for cluster X-ray luminosity evolution, and 
Mullis et al. (2001) described the NEP supercluster discovered 
in the survey X-ray data. The extragalactic component of the 
NEP survey was discussed in Mullis et al. (2003, 2004a,b), 
Gioia et al. (2004), and Henry et al. (2006). 

The identification of the optical counterparts (based both 
on optical photometry and spectroscopy of all candidate coun- 
terparts) to each of the X-ray sources was completed by Gioia 
et al. (2003). This X-ray based data set is unique in its com- 
bination of contiguous sky coverage (80 deg^), and sensitivity, 
which is higher than any of the Einstein-hased surveys. 

More recent X-ray surveys, such as the Chandra Deep 
Fields, are far deeper than the NEP survey, but have small sky 
coverage. These surveys are biased towards older stars, which 
are more numerous but fainter in X-rays, with the younger stel- 
lar population being very poorly sampled. The small sky cov- 
erage does not, however, allow robust conclusions to be made. 

The stellar content of the Chandra deep survey (Feigelson 
et al. 2004) reveals that the old star population does not contain 
any excess of yellow stars, but instead a deficiency of dG-dK 
stars, while the signatures of dM stars are reproduced perfectly 
by galactic models. This behaviour is the opposite of that ob- 
served for the young stellar population, and suggests that some 
changes occur at intermediate ages, where NEP survey may 
give a unique contribution thanks to its intermediate sensitiv- 
ity. Thus the NEP survey is unique in its ability to properly 
sample the young and intermediate-age stellar populations in 
the Galactic disk. 

In Paper: we analyzed the stellar content of the NEP sur- 
vey, to determine the nature of active stars in the solar neigh- 
borhood. We determined the spectral types for all source detec- 
tions. We then compared the observational data with the pre- 
dictions of the XCOUNT model (Favata et al. 1992; Micela 
et al. 1993), and, at intermediate X-ray fluxes, ascertained the 
presence of an excess of active yellow stars, in the solar neigh- 
borhood. 

To characterize the stellar population of the NEP survey 
and determine the nature of the observed excess, we derived, 
in the second paper of this series, the physical properties of 
the X-ray sources, such as radial and rotational velocities, and 
lithium abundances. In a forthcoming paper, we will determine 
the chemical composition of the stellar photospheres; this will 
aUow the study of the enrichment history of the interstellar mat- 



Table 1. Sampled stars per spectral type 



Sp.Type 


N.Obs.(Paperi) 


N.Obs.(This work) 


A 


3 


1 


F0-F5 


10 


9 


F6-F9 


8 


4 


G 


29 


18 


K 


53 


22 


M 


41 


2 


Total 


144 


56 



ter and, combined with kinematical data and ages, will provide 
a powerful way of probing the chemical and dynamical evolu- 
tion of the Galaxy. 

The present paper is organized as follows: in Sects. 2 and 3, 
we describe the selection rules of our sample, the observations, 
and the data reduction. In Sect. 4, we describe the methods of 
analysis and discuss the distribution of radial and rotational ve- 
locities, while binary stars and lithium-abundance determina- 
tion are described in Sect. 5 and 6. In Sects. 7 and 8, we discuss 
and summarize our findings. 

2. Observations and data reduction 

To complete our proposed investigation, we would have ide- 
ally wanted to observe the entire sample of NEP stars. Due 
to telescope and instrumentation constraints we instead had 
to restrict our sample to stars with J < 11 in the 2MASS 
Catalogue; this selection criterion should not affect the achieve- 
ment of our goals, apart from the study of dM stars. In Table[T] 
we summarize the spectral type classification of stellar X-ray 
NEP sources, as derived in Paper i, together with the sample 
observed in the present work. As can be seen from Table [1] 
we covered almost the entire sample of early-F (F0-F5) stars, 
while we sampled a fraction between 40% and 60% of late-F 
to G and K stars; dM stars were poorly sampled (due to their 
faintness at the energies considered) and we are unable to mea- 
sure reliably their mean properties, even though, with dK stars, 
they form a substantial part of the data set. Finally we can state 
that our sample is representative of stars earlier than dM. 

The observations, with signal-to-noise ratios between 40 
and 440, were obtained (in service mode) during several ob- 
serving runs between January and June 2004 (see Table lA.U 
at the Telescopio Nazionale Galileo (TNG) located at the 
Observatorio del Roque de Los Muchachos (La Palma, Canary 
Islands). 

Our survey contains 56 stars in the region of sky of 80.6 
deg2 about the North Ecliptic Pole (NEP) at q'=18''00'"00', 
5=-h66°33'00" (J2000.0), at a Galactic latitude, b = 29.°8, with 
spectral types in the range F0V-M3V, and values of B- V in the 
range 0.3 to 1.5, which corresponds approximately to effective 
temperatures of about 7300 K-3500 K. For all observations, we 
used the cross-dispersed echelle SARG spectrograph attached 
to the 3.58 m TNG. The dioptric camera has a focal ratio of 
//4.8, and the detector is a mosaic of 2 2k x 4k thinned and 
back illuminated CCDs, with a 13.5 /zm pixel size, and a scale 
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of 0.16 arcsec/pixel. We used the Yellow grism, with a disper- 
sion of 61.0 A/mm, the FW4 filter, and selected a binning of 
2x1. The spectra continuously cover a wavelength range from 
4620 A to 7920 A, with a resolving power of about 57000, 
yielding a FWHM of 5 km s ^ 

To enable a differential analysis, we obtained a solar flux 
spectrum, i.e. a high S/N 500) moonlight spectrum acquired 
using the same instrumental configuration as the target stars. 

As preparation for our cross-correlation analysis, the objec- 
tive of the data reduction is to produce wavelength-calibrated, 
one-dimensional spectra for each echelle order, free of relevant 
instrumental and telluric signatures. Tasks, implemented in the 
IRAlQ package, were used to carry out the standard reduc- 
tion procedures on echelle spectra (i.e. bias subtraction, flat- 
field correction, removal of scattered light, order extraction, 
and wavelength calibration). The two CCDs were separated 
and reduced independently. Numerous emission spikes, possi- 
bly arising from CCD electronics, were initially found to affect 
the extracted spectra. To remedy this problem, the data were 
extracted after excluding deviant pixels in the two-dimensional 
science frames. The rejection thresholds of the sigma-clipping 
algorithm were varied to remove the artefacts efficiently, with- 
out affecting significantly the remainder of the spectrum. The 
wavelength calibration was performed using Thorium-Argon 
lamp exposures acquired each night. For three nights, Thorium- 
Argon lamp exposures were not acquired; for data obtained 
during this time, wavelength calibration was completed using 
lamp exposures available for the night closest in time; we iden- 
tify the spectra affected by labelling the corresponding stars 
with an asterisk in Table lA.fl (we use a similar symbol in all ta- 
bles). We compared the wavelength shifts between telluric lines 
of these spectra, and of spectra for which wavelength calibra- 
tion was completed using lamp data acquired during the same 
night; we corrected each spectrum by the corresponding wave- 
length shift (which amounts, at most, to 0.9 km s '). 

3. Analysis methods and results: radial and 
rotational velocities 

In Paper i, we analyzed the stellar content of the NEP survey 
and we determined spectral types of all of the stellar coun- 
terparts to the 144 ROSAT NEP X-ray sources. Our analysis 
confirmed the presence of an excess of active yellow stars, in 
agreement with previous results for shallower surveys. To char- 
acterize the nature of the NEP stellar sample, we determine the 
radial velocity (RV) of each star, which provides a kinematic 
index of the population to which they belong, and the projected 
rotational velocity vsini, which is strongly correlated with the 
X-ray luminosity in single stars (Pallavicini et al. 1981), and 
evolves with stellar age (Skumanich 1972). 

We derive simultaneously the radial and rotational veloc- 
ities using the Fourier Cross-correlation technique originally 
described by Tonry & Davis (1979). The cross-correlation tech- 
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' IRAF (Image Reduction and Analysis Facility) is distributed 
by the National Optical Astronomy Observatories, operated by the 
Association of Universities for Research in Astronomy, Inc., under 
cooperative agreement with the National Science Foundation. 



Spectral Type 

Figure 1. Heliocentric radial velocity distribution of the NEP 
sample stars with respect to spectral type (temperature). We 
assumed that the weighted mean of three independent measures 
is our estimate of radial velocity for each star. The uncertainty 
was computed by combining in quadrature the uncertainties of 
the measures. 

nique was developed to derive robust radial velocities using all 
of the information in the spectrum at once, rather than one spec- 
tral line at a time. We coiTelate the spectrum of a program star 
against that of a stellar template. The position of the correla- 
tion peak indicates the relative radial velocity, and the breadth 
of the peak is interpreted as being due to rotational broadening. 

3.1. Radial velocities 

To check the robustness of our approach, we compared the 
cross-correlation method with that performed line by line, then 
we calculated radial velocities adopting three different ap- 
proaches which rely on: 

1 . Wavelength shifts of selected photospheric absorption fea- 
tures 

2. Cross-coiTelation function of the object spectrum with that 
of a template star belonging to the sample 

3. Cross-coiTelation function of the object spectrum with that 
of the Sun, obtained with the same instrumental configura- 
tion as the target stars 

We comment on the application of all three methods: 

1 . The heliocentric RVs of the sample stars were obtained by 
measuring the wavelength shifts of a set of 58 measurable 
and unblended prominent photospheric absorption features 
(Table 111, using the IRAF task RVIDLINES, and adopting 
the weighted mean of the results (using the inverse of the 
variances as weights). The variance of the weighted mean 
is the inverse of the sum of the weights. 

2. We calculated the RV of the sample stars using the IRAF 
task FXCOR, which performs a Fourier Cross-Correlation 
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Table 2. Absorption features selected to calculate radial veloc- 
ities from wavelength shifts 

Wavelength (A) Element Wavelength (A) Element 

6419.95 Fei 6705.10 Eel 

6421.35 Fei 6710.32 Eei 

6430.85 Eei 6717.69 Cai 

6439.07 Cai 6726.67 Eei 

6449.81 Cai 6743.12 Tii 

6450.23 Co I 6750.15 Eei 

6455.60 Cai 6767.77 Nil 

6462.57 Cai 6783.70 Eei 

6469.19 Eei 6786.86 Eei 

6471.66 Cai 6806.84 Eei 

6475.62 Eei 6810.26 Eei 

6481.87 Eei 6814.95 Coi 

6482.80 Nil 6820.37 Eei 

6493.78 Cai 6828.59 Eei 

6494.98 Eei 6839.83 Eei 

6495.74 Eei 6841.34 Eei 

6496.47 Eei 6842.04 Nil 

6496.90 Ban 6842.69 Eei 

6569.22 Eei 6843.66 Eei 
6572.78 Cai 7090.38 Eei 

6574.23 Eei 7095.43 Eei 
6575.02 Eei 7110.90 Nil 
6581.21 Eei 7112.18 Eei 
6586.31 Nil 7122.20 Nil 

6592.91 Eei 7130.92 Eei 
6593.87 Eei 7142.52 Eei 
6598.60 Nil 7145.32 Eei 
6609.11 Eei 7148.15 Cai 
6703.57 Eei 7151.49 Eei 



of a given spectrum with a "template" stellar spectrum, 
as described by Tonry & Davis (1979). For each star ob- 
served, we cross-correlated its spectrum with that of the 
template using the wavelength range between about 6625 
A and 6710 A, which contains no significant telluric ab- 
sorption lines. While performing the cross-coiTelations for 
individual stars, the wavelength ranges used were visually 
inspected and obvious emission lines, cosmic rays, and/or 
bad pixels were masked. As our stellar template, we used 
a high signal-to-noise (S/N x 180) spectrum of the bright 
star nep 6163 (V ^ 7.6), which has a modest projected rota- 
tional velocity (vsini ^ 3.7 km s"', derived by comparison 
with synthetic spectra computed using ATLAS9). The he- 
liocentric RV of the template object nep 6163 was obtained 
by measuring the wavelength shifts of a list of prominent, 
photospheric, absorption features, as for method 1 . The he- 
liocentric RV of each observed object is then calculated by 
adding the heliocentric RV of the template star, to the RV of 
each observed star, with respect to the template. The uncer- 
tainty in each value of heliocentric velocity was computed 
by combining (in quadrature) the uncertainties in the helio- 
centric RV of the template star, and in the relative RV of 
the object, measured by cross-correlation with the template 
star. 
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Figure 2. Histogram of the heliocentric radial velocities, the 
Sun velocity with respect to the Local Standard of Rest (LSR) 
is indicated (\visr\-13.4 km s"', Binney & Merrifield, 1998). 
The superimposed gaussian radial velocity distribution for the 
targets with a standard deviation of 13 km s"', is shown. 



3. We performed an analysis similar to method 2, using as a 
stellar template, a high signal-to-noise (S/N a: 400) solar 
spectrum, acquired with the same instrumental configura- 
tion as the sample stars. Using IRAF task RVCORRECT, 
we find the heliocentric coiTection to the observed radial 
velocity of the solar spectrum, due to the relative motion 
Sun-Earth. To measure the RV of each observed star, we 
cross-correlated each spectrum, with respect to the cor- 
rected template spectrum. 

The three measurements agree to within a few km s the 
maximum difference being about 2.5 km s we therefore use 
the weighted mean, which is obtained using the inverse of the 
variances as weights, as our radial velocity measurement for 
each star; the variance that we associate with this weighted 
mean is the inverse of the sum of the weights. 

For a few stars observed, the peak of the cross-correlation 
function is asymmetric, probably due to binarity and/or high 
rotation rates (> 60 km s"'), particularly in early F-type stars; 
we are unable however to determine whether the stars are 
binaries and/or fast rotators because our spectra were obtained 
during a single night. For about 7 stars, this asymmetry is 
indicative of a binary nature. These stars are discussed later 
Figure [T] shows the heliocentric radial velocities as a function 
of spectral type for the whole sample. The histogram in Fig.|2] 
provides the radial velocity distribution for single stars, which 
peaks significantly, with most velocities being in the range 
between and -30 km s"', and has a mean value of about -15 
km s"^ 
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Figure 3. Calibration curve used for the vsini estimation from 
the width of the cross-correlation peak. The curve was created 
by cross-correlating the spectrum of the template star with a set 
of synthetic spectra with different rotational velocities. 



Figure 4. Histogram of the rotational velocities; almost all of 
the sample stars have rotational velocities that are comparable 
or higher than that of the Sun (x2 km s '). 



3.2. Rotational velocities 

The cross-correlation function measures the projected rota- 
tional velocities (vsini) of our targets, relative to stellar spectra 
of essentially zero rotational velocity (our "templates"); this is 
possible because the measured width of the cross-correlation 
peak is a function of the broadening of the spectrum lines, 
due to rotation velocity. Assuming that the line broadening of 
the NEP stellar spectra is dominated by rotational broadening, 
and that our stellar template is an intrinsically narrow-lined 
spectrum, the width of the cross-correlation peak is a mea- 
sure of the vsini of the NEP targets, which may be calibrated 
as discussed below. We use a standard Local Thermodynamic 
Equilibrium (LTE) spectral analysis with the latest version of 
the spectral-line analysis code MOOG (Sneden 1973), and a 
grid of ATLAS9 model atmospheres (Kurucz 1993), to com- 
pute synthetic spectra for comparison with the observed spec- 
tra. 

The ATLAS9 stellar atmospheres are computed without the 
overshooting option, and by assuming a mixing length to pres- 
sure scale height ratio of a=0.5. Assumptions made in the mod- 
els include: the atmosphere is plane-parallel and in hydrostatic 
equilibrium, the total flux is constant, the source function is de- 
scribed by the Planck function, and the populations of different 
excitation levels and ionization stages are governed by LTE. 
The atomic data used to generate synthetic spectra are taken 
from Kurucz Atomic Line Database (Kurucz 1995). 

The cross-correlation is performed with FXCOR. The main 
peak of the cross-correlation function is fit using a Gaussian 
function and its FWHM is measured. These FWHM values 
are converted to vsini values, using the following calibra- 
tion process. We compute a grid of synthetic spectra, assum- 
ing logg=4.5 (cm s"^), microturbulent velocity ^=1.0 km s"' 



and solar metallicity, which are reasonable parameters for our 
template star, assuming the SARG instrumental resolution, 
R=57000. We convolved these spectra with a rotational broad- 
ening function with vsini varying from to 50 km s"'. We al- 
lowed Teff to vary between 4000 K and 7500 K (which are the 
typical temperature values for our targets) in 500 K intervals, 
to check the sensitivity of this relationship to temperature. 

The synthetic spectra were cross-correlated with our tem- 
plate and the FWHM of the cross-correlation peak was mea- 
sured. We compared the results derived for the same analysis 
using synthetic spectra computed for \ogg-2>.Q (cm s"^). The 
comparison highlights that the calibration results are almost in- 
dependent of gravity, while the dependence on temperature is 
more pronounced. 

The resulting relation between the FWHM of the cross- 
correlation peak and the corresponding vsini, for the lowest 
(Feff = 4000 K) and the highest {T^s = 7500 K) temperature 
values, are shown in Fig. [3] 

We note that the discrepancies between the calibration 
curves, calculated for different T^s, are consistent within the 
typical errors inferred for the rotational velocities measured 
for our targets. Using the relationship corresponding to the ex- 
pected temperature of our template (according to spectral type), 
we linearly interpolated the FWHM of the cross-correlation 
peak, with the derived calibration curve to estimate the vsini 
values. 

The errors of the rotational velocities were calculated as pre- 
scribed by Rhode, Herbst & Mathieu (2001), to be vsini I {I +r), 
where r is the Tonry & Davis (1979) parameter, which is a mea- 
sure of the signal-to-noise ratio of the cross-correlation peak. 
When the calibration derived is applied to determine vsini for 
each star, we must pay attention to the minimum value of vsini 
to which this procedure is sensitive. The measured width of the 
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cross-correlation function can be described by the following 
formula (Tonry & Davis, 1979): 

2 -2,2,^2 

where cr„,f is the rotational width, cr„a, is the natural, or intrin- 
sic, line width which is negligible, and cr,>,j, is the instrumental 
broadening (the last term of the equation includes two contribu- 
tions because two spectra form the cross-correlation function). 

For the derivation of the minimum detectable vsini, we fol- 
low Bailer-Jones (2004). The terms are additive; with perfect 
data of "infinite" S/N ratio, we could therefore determine "any" 
rotational broadening, even with non-zero instrumental broad- 
ening. In practice, this is not however possible due to noise. 
Bailer- Jones (2004) assumes that cTroi is detectable only if its 
value exceeds that of other broadening contributions, and its 
minimum measurable value is V2cr,„„. Our theoretical FWHM 
of the instrumental broadening is fixed by the slit width to be 5 
km s We note that the full width of a rotational profile cor- 
responds to twice the rotational velocity: one half of the line 
is created by the blueshifted approaching limb of the star, the 
other half by the redshifted receding limb (Bailer-Jones, 2004). 
The minimum detectable vsini for our data is V2 x 5.0/2 = 3.5 
km s ' , which corresponds to the flat portion of the calibrations 
in Fig. [3] 

The cross-correlation technique measures the projected, ro- 
tational velocity of the Sun to be 1.63+0.34 km s"'. However, 
this measurement, which agrees with the well-known value 
assumed for the Sun (^ 1.86 km s"', Soderblom 1982), is 
below the minimum detectable vsini of our study. As we can 
see from Fig. [3] this value appears in the flat portion of the 
calibration relation; the adoption of the mean instrumental 
broadening as the upper limit on the vsini for the slowest 
rotators, is therefore conservative. 

We indicate the vsini values derived for slow rotators with < 
3.5 km s"'. 

To perform a further comparison, we adopted a second 
method to calculate rotational velocities. 
We performed a detailed comparison between the observed 
and synthetic spectra in the wavelength range 6420 A - 
6480 A, which includes several Fe I and Ca I absorption 
lines, of moderate strength (» 100 mA). Using Kurucz model 
atmospheres and the MOOG software, we generated a grid 
of LTE synthetic spectra. We allowed the stellar temperature 
vary, and assume log^=4.5 (cm s~^), ^=1.0 km s"', a solar 
metallicity, and a macroturbulent velocity, that is constrained 
by the spectral type (Fekel, 1997); the rotational velocity 
was then the only free parameter. Each spectrum was finally 
convolved by a 0. 1 A/px Gaussian FWHM, which reproduces 
the instrumental profile, as estimated by our lamp spectra. For 
the determination of rotational velocities, we used the neutral 
iron fines at 6430.85 A, 6546.24 A and 6569.22 A and the 
Cai lines at 6439.07 A and 6471.66 A, which were strong 
but not saturated lines, and unblended with other lines. We 
compared the observed and synthetic spectra and determined 
by eye the value of vsini that provides the most appropriate 
fit to the observed target rotational velocity. We assumed that 



the spread in vsini values for individual lines represented 
the uncertainty in each measurement. We note that the vsini 
values were measured by comparing observed spectra and 
synthetic profiles that had been convolved by all broadening 
effects. The minimum value of vsini to which this procedure 
is sensitive is therefore formally zero. The two independent 
measurements agree to within a few km s"'; the maximum 
diff'erence is about 3.1 km s"' (for the low S/N spectra). We 
therefore calculated the weighted mean of the two measure- 
ments of vsini and the inverse of the sum of the weights 
was taken to be the variance associated with the weighted 
mean. Figure|4]shows the histogram of the rotational velocities. 



4. Binary stars 

For approximately 7 stars, ^ 12% of all stars observed, the 
shape of the FXCOR cross-correlation function peak was 
asymmetric; this may be indicative of a binary nature for these 
stars. The peak of the cross-correlation function can be fitted 
using two overlapping Gaussians, instead of one single func- 
tion, when the blending between the two stars is not severe; 
an example of such a case can be seen in Fig. |5] However, we 
cannot measure the radial velocity of binary systems, because 
our spectra were acquired during a single night. Both radial 
and rotational velocities for such stars, when measurable, refer 
to the contributions of both stars. We report these results for 
completeness and stress that they are affected by blending of 
both contributions, and must therefore be taken with extreme 
caution. The only exception is the spectroscopic binary (SB2) 
nep 5520, for which spectral lines from both components 
are clearly separated and the calculations refer separately 
to each star. For this system, the two components appear to 
have similar spectral types and eff'ective temperatures, i.e. to 
contribute equally to the continuum. 
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Figures. Example of peaks of the cross-correlation function 
produced by cross-correlating a binary with our template star. 
Relative radial velocities of the two components are indicated 
in the panel. 

5. Lithium abundances 

Of the 56 total observed stars, 25 show a prominent lithium 
feature, 5 of which are supposed to be binaries. For these stars. 



nepaaao - ko 




L. Affer et al.: The stellar North Ecliptic Pole Survey 

Table 3. Properties of high radial velocity stars 



6710 
Wavelength (A) 



Figure 6. Spectra in the lithium region for two NEP sources: 
nep 3220 does not show the lithium line while nep 4780 shows 
a prominent line. The spectra have been continuum subtracted 
and offset on the vertical axis. 

we used the IRAF task Continuum to perform the spectrum 
normalization. We created a synthetic model atmosphere, us- 
ing the ATLAS9 code (Kurucz 1993), assuming \ogg-A.5 (cm 
s"^), ^=1.0 km s"' and solar metallicity, and derived the effec- 
tive temperatures using the star's spectral type. 

These models were used to create synthetic spectra for a 
small interval of approximately 50 A about the lithium feature 
at 6707.8 A, from 6700 A to 6750 A, to determine the Hne-free 
regions that were fitted using low-order cubic spline functions, 
and a variable residual rejection limit. Both the line-free spec- 
tral range and the rejection limit were chosen interactively in 
each case, based on the visual inspection of fitting results. In 
Fig. |6] we show the normalized spectra about the Lii 6707.8 
A line for two stars observed: the spectrum of nep 3220 does 
not show the lithium line while nep 4780 shows a prominent 
line. The Li i 6707.8 A fine is close to the Fe i 6707.44 A line 
and both lines are often blended, due to spectral line broad- 
ening by stellar rotation. The equivalent width of each of the 
partially-blended lines, Li i 6707.8 A and Fe i 6707.44 A, was 
determined using a two-Gaussian fitting routine provided by 
IRAF. To estimate uncertainties in the EW measurements, we 
repeated both the continuum normalization and EW measure- 
ment three times. We then accepted the weighted mean of the 
three measures to be final EW measurements and as typical er- 
ror the standard deviation of eiTors, which were typically not 
more than a few percent (^3 mA). 

For stars with 4000 K <reff<6500 K, the Li abundances 
N(Li) were derived using the measured EWs and from a lin- 
ear interpolation of the growth curves of Soderblom et al. 
(1993). Effective temperatures were derived from the spectral 
type (Allen, 1973), as determined in Paper:; we derived the er- 
ror in two subtypes from the published spectral types (Paper:), 
and this uncertainty was propagated into the Li abundances. 



nep 


Sp.Type 


High Vrad 


High Vrot 


Li 






(> |30| kms-') 


(> 15 kms-') 




1690 


G2 


-37.0+2.7 






1800 


F5 


-58.0±4.1 


16.3+1.4 


3.30+0.20 


3560 


G8 


-26.0±1.3 




2.30+0.30 


*4090 


G2 


-41.6+2.2 






4810 


F5 


-37.0+1.6 






5510 


KO 


-53.5+3.0 


20.4+2.1 


3.70+0.50 



*: See NOTE in Table EH 



The lithium EW for binary stars was computed in rela- 
tion to the combined continuum, which was the sum of con- 
tinuum contributions from both components. Without informa- 
tion about the individual spectral types of both components, we 
were unable to separate the continuum contributions and obtain 
individual estimates of the lithium EW of each star (apart from 
the SB2 nep 5520). We measured the EWs for the suspected 
binary stars, but we emphasize that these measurements corre- 
spond to the unsolved blend including the contribution of both 
stars; we can therefore only state that these stars have lithium. 
While this is the best one can do with the available information, 
this procedure is likely to introduce additional eiTors. Therefore 
in Fig.|7] in which we report the lithium EWs measured for the 
subsample of single stars, such system are flagged with a dif- 
ferent symbol (filled squares, while the EWs for the two com- 
ponents of nep 5520 are indicated as ' x'). 

In Table lA.2l we indicate the Li abundances N(Li) derived, 
together with the NLTE corrections interpolated from the grid 
of Carlsson et al. (1994). When no lithium line was visible in 
the spectrum, we used the statistical noise measured from the 
spectral continuum, to determine an upper limit to the Li line 
EW, except for those stars with such high rotational velocity, 
that the fines were completely smeared out. To set constraints 
on the age of our targets, the lithium EWs of the Pleiades (10^ 
yrs) and Hyades (10^ yrs) clusters are plotted in Fig.|7]for com- 
parison (Soderblom et al., 1993 and 1990, respectively). 

The comparison of the three samples of stars shows that 
the EWs of most targets lie between those of stars of similar 
temperature in the two clusters, and in a few cases are lower 
than those of the Hyades. 

6. Discussion 

We derive both radial and rotational velocities, and lithium 
abundances, for the stellar population of the NEP survey. 

We use different methods to determine the same parame- 
ters, which enables us to infer the reliability of the analysis 
techniques employed. We obtain results that are consistent with 
each other, which demonstrated the robustness of our analysis. 
We quote the results in Table IA.2I col. 1 identifies stars with 
NEP numbers; col.s 2-3 provide the heliocentric RV and the 
vsini values, with estimated eiTors (we indicate that 3.5 km s"' 
is the upper limit for measurements of rotational velocity that 
are below the instrumental profile width, as discussed in Sect. 
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|3.2| i; cols. 4-7 indicate the equivalent widths of the lithium line 
(6707.8 A), the temperature and its associated error (derived 
from spectral type), and lithium abundances (derived from the 
interpolation of the growth curves of Soderblom et al., 1993), 
respectively; col. 8 reports the spectral types and in col. 9 we in- 
dicate whether the star is a binary (label "x") or a double-line 
spectroscopic binary (label "SB2"). Label "*" indicates that 
we lack the Thorium-Argon lamp exposures for the nights in 
which data for these stars were acquired. These observations 
were wavelength-calibrated using lamp exposures acquired for 
the night closest in time. We assumed the same spectral type 
and temperature for the two component of SB2 nep 5520; this 
is an "ad hoc" assumption, which is reasonable because the 
lines are similar In Table lA.2l we indicate our poor knowledge 
about binaries with the label "?", close to the spectral type in- 
dication. 
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Figure 7. Distribution of lithium equivalent widths of single 
NEP stars (filled dots) and binary stars (ringed dots; 'x' for 
nep 5520), with respect to spectral type, compared to Pleiades 
(empty squares) and Hyades (empty dots) distributions. Arrows 
indicate EW upper limits for those stars in which no lithium 
line was visible in the spectrum. 



6.1. Radial velocity 

In Figure |2] we show the histogram of the heUocentric radial 
velocities, and superimpose a Gaussian radial velocity distri- 
bution for the target stars, adopting a standard deviation of 13 
km s that corresponds to the histogram variance, and a mean 
value of -13.5 km s"', that agrees well with the bulk of the data. 

The Sun moves relative to the Local Standard of Rest 
(LSR) with a mean velocity of |vl5r|=13.4 km s"' (Binney & 
Merrifield, 1998). 

The Gaussian radial velocity distribution, which represents 
a good fit for most of the targets reflects the motion of the Sun 




Figure 8. Rotational velocity distribution with respect to spec- 
tral type (temperature); average trends of rotational velocity 
of Pleiades, Hyades and field stars, are shown (Favata et al. 
1995; Kraft 1967; Soderblom 1983a,b; Skumanich 1972). With 
squares we indicate stars that have Nai(Z)i,D2) emission; "x" 
indicates the presence of Ha emission and diamonds are for 
stars that have both Ha and Nai(Di,D2) emission. 



relative to the LSR. The distribution of radial velocity is consis- 
tent with that of a young population of disk stars of age 4x10^ 
yrs (Wilson 1963; Kraft 1967; Skumanich 1972; Simon 1985), 
or younger. We detect a small number (6) of stars with large ra- 
dial velocity (> |30| km s"^). Nevertheless we cannot consider 
these stars as a statistical meaningful sample, because only 4 
of these stars are within 2cr from the mean value, one is within 
3cr and only one has a velocity which falls beyond 3cr from the 
mean value with an extremely small probability (|P| > 3 cr x 
0.3%). 

In Table [3] we report the properties of these stars and the 
lithium abundances. In particular we want to know if these stars 
are older than the remainder of the observed stars or if they are 
peculiar in some way. Of these six stars, three of them (nep 
1690, nep 4090, nep 4810) are slow rotators (< 15 km s"') that 
have no detected 6708 A lithium line. These could be old stars. 
One of these stars (nep 3560) is a slow rotator and a Li line is 
detectable in its spectrum, and finally, two stars (nep 1800, nep 
5510) are fast rotators (> 15 km s"') and have measurable Li 
line in their spectra. 

A possible explanation for the "high velocity" stars can be 
given in terms of binarity. Nep 5510 is probably a binary sys- 
tem. Analyses of spectra acquired for all other stars do not, 
however, confirm conclusively the presence of binarity. 

Abundance analyses, in particular the measurement of the 
metallicity and chemical abundance ratios of several elements 
([a/Fe], [O/Fe]), could provide insight into the stellar popula- 
tions, because abundances are a key indicator of which popula- 
tion a star belongs to. 
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6.2. Rotational velocity 

Figure |4] shows that almost all stars have projected rotational 
velocities, which are interpreted as lower limits to actual val- 
ues, that are comparable to or higher than that of the Sun 2 
km s"'). 

Rotational velocity decreases with increasing age. Since 
magnetic activity depends on rotation and differential rotation, 
the level of magnetic activity and hence X-ray emission will be 
different for fast and slow rotators. 

Our sample seems to be dominated by relatively young or 
intermediate-age stars. This statement is supported by Fig. [8] 
which shows the vsini distribution as a function of the spectral 
type and compares it to the average trends for the Pleiades, 
Hyades and field stars; we indicate stars that have Nai(Di,D2) 
emission with square symbols. Ha emission with x symbols, 
and stars that have both Ha and Nai(Di,D2) emission with 
diamond symbols. We see that most stars have measurements 
that are consistent with an age between 10** yrs (Pleiades age) 
and 4x10^ yrs (field stars), compared to the age of the Sun, 
5x10'' yrs. 

A small number of K stars have measured rotational veloc- 
ities that are higher than measurements for the Pleiades. We 
propose that these stars have ages that are comparable to those 
of the Pleiades, or are even younger. Two of the stars could be 
binaries. For the remainder of the sample the youthness is the 
most plausible explanation. 

Widely-separated binaries behave in a similar way to sin- 
gle stars because tidal forces can be neglected (Slettebak 1963, 
Weiss 1974). Close binaries, however, show altered rotation. 
Close binaries, of early spectral type, rotate more slowly than 
single stars, while those of later type rotate much faster than 
their single counterparts (Gray 1976, 1984). 

In Figure s|9] and [To] we report the LTE lithium abundances, 
upper limits included, as a function of rotational velocity and 
temperature, respectively. Stars that show high lithium abun- 
dances and low rotation are probably young; this hypothesis is 
reinforced by the presence of emission in the Ha and Nai Di 
and D2 lines, in most spectra (see Table lA.3l l. 

6.3. Comparison with model results 

The predictions of the model XCOUNT (Favata et al. 1992; 
Micela et al. 1993) for the complete NEP survey, for each spec- 
tral type and age range are summarized in Table |4] 

XCOUNT assumes an exponential disk-like spatial distri- 
bution of the stars, with a radial scale length of 3.5 kpc, as in 
Bahcall & Soneira (1980). We followed the modifications of 
the Bahcall & Soneira (1980) model of Micela et al. (1993), 
which introduced an age-dependent scale height, dividing stars 
into three age ranges: 10^-10^ lO^-lO'', and lO'^-lO'" years, 
with scale heights of 100, 200 and 400 pc, respectively. Micela 
et al. (1993) assumed X-ray luminosity functions derived from 
ROSAT observations of the Pleiades (Micela et al. 1996) and 
Hyades (Stern et al. 1995), and from Einstein data for nearby 
stars (Schmitt et al. 1985; Maggio et al. 1987; Barbera et al. 
1993), which are considered to be prototypes for the three age 
ranges. 




Figure 9. Lithium abundances as a function of the rotational 
velocity: diamonds indicate stars with vsini greater than those 
of Pleiades; "x" indicate vsini values between the average 
curves of Pleiades and Hyades; dots indicate vsini values be- 
tween the average curves of Hyades and field stars; down- 
pointing arrows indicate upper limits of lithium abundances. 
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Figure 10. Lithium abundances as a function of the effective 
temperatures: symbols as in Fig. |9] We indicate Pleiades and 
Hyades fits, obtained by Favata et al. (1993). 



We scaled the XCOUNT predictions by a factor corre- 
sponding to the fraction of the stars that we observed (^ 0.6 
for late-F to G stars and ^ 0.4 for K stars), and assigned each 
star to a given age range based on its lithium EWs and/or rota- 
tional velocities. We used lithium as a primary indicator and as- 
signed a star: to the 10^-10** years age range, if its lithium EWs 
were greater than (or comparable to) those of the Pleiades; to 
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Table 4. Summary of X-ray source counts predictions for each 
spectral type and range of age derived from XCOUNT (from 
Paper:). 



Sp.Type 


N(young) 


N(intermed.) 


N(old) 


Total 


A 


0.63 


1.41 


0.28 


2.32 


F0-F5 


0.32 


1.53 


7.34 


9.19 


F6-F9+G 


1.91 


8.42 


7.63 


17.96 


K 


4.24 


8.86 


20.26 


33.36 




7.1 


20.22 


35.51 


62.83 



the 10^-10^ years age range, if its lithium EWs were between 
the measurements shown by the average curves of Pleiades and 
Hyades; and to the 10^-10'" years age range, if its lithium EWs 
were less than those of the Hyades. For those stars for which 
lithium is not detected, we used rotational velocities as an al- 
ternative indicator, applying the same criteria as for the lithium 
EWs. We focus on late-F to G and K type stars because our age 
indicators (lithium and rotation) are not appropriate for early- 
type stars. For late-F to G stars, we expect to find about 1 .2 stars 
that have an age < 10*^ yrs, but instead find 7 stars. Similarly, 
we expect about 5 stars to have an intermediate value of age 
in the range 10^-10^ years, and indeed we find 4 stars. For the 
age range 10^-10'" years, we find 6 stars in contrast to about 5 
predicted by the model. For K stars, we find 7 stars that have 
an age < 10^ yrs, compared with about 2 stars predicted by the 
model; 5 stars to have an intermediate age, between 10**- 10'' 
years, with respect to about 4 predicted stars; and finally 7 stars 
in the range between lO^-lO'" years, in contrast to a prediction 
of 8 stars. Our results are summarized in Table |5) we indicate 
the observed and predicted (by XCOUNT model) stars for our 
sample, and, between brackets, the total counts of stars that 
may have an age in the range considered, on the basis of ro- 
tational velocities and lithium EWs, including those stars with 
only rotational velocity and no Uthium measurements. For 8 of 
the late-F to G and K stars observed, we did not measure the 
upper limits of lithium EWs due to the low S/N ratio of the 
spectra and/or high rotation. 

The observed excess seems to be concentrated at young 
age. Our comparison suggests a burst of star formation, with 
an increasing of new-born stars by a factor of » 4.5 (or « 3.5 if 
we exclude 4 probable binary systems, concentrated at young 
age) in the last 10^ years, while it is consistent with a constant 
stellar birthrate in the intermediate and old ages, between 10'^- 
lO'" years. 

7. Conclusions 

We derived fundamental quantities, such as radial and rota- 
tional velocities and lithium abundances, for stars detected in 
the NEP X-ray survey. In Paper i, we confirmed the presence of 
an excess of yellow stars in the NEP sample and suggested that 
such an excess could be due to the presence of a young popu- 
lation, formed by a recent event of star formation. Moreover, 
there was no evidence for an equivalent excess of M stars, 
which would be necessary to identify the population as really 



young. We refer the reader to Paper: for a possible explana- 
tion of the absence of an excess. In the present work, M stars 
are poorly sampled (^ 0.05% of the entire sample of NEP M 
stars) and, as a consequence, we cannot state anything about 
their mean properties. Our results, which refer to F, G and K 
stars, confirm that most of the stars seem to belong to a young 
age population. The population excess of the model prediction 
described in Paper: appears to be composed of mainly young 
stars (with ages comparable to those of Pleiades). 

Our data suggest a burst in the birthrate, by a factor of 
X 4, in the last 10*^ years. The Gould Belt is an asymmetry 
in the apparent distribution of the brightest stars in the sky, 
with respect to Galactic equator, probably associated with 
a recent star formation episode close to the Sun (Guillout 
et al 1998). The NEP surveyed area is far from the Gould 
Belt and this star formation episode was not accounted for 
in the XCOUNT model, and is possible that other local, 
less prominent events of star formation occutTed in the solar 
neighborhood. Nevertheless few stars, among the fastest K 
rotators, may be old binary systems with tidally locked rotation. 
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Appendix A: Observing log for the NEP spectra 
and results 

In this section, we present the observing log for the observa- 
tions of stars detected in NEP survey (hereafter NEP spectra) 
and the results derived during the present study. 
In Table lA.ll we report the observing log for the NEP spec- 
tra, including identification numbers, observation dates, right 
ascension, declination (Eq. 2000.0), spectral types (as derived 
from Micela et al. 2007), and total observing times. 
In Table IA!21 we provide our measurements: heliocentric RV 
and vsini values, with estimated errors, equivalent widths of 
lithium line (6707.8 A), temperatures with associated errors 
(derived from spectral type), and lithium abundances (derived 
from the interpolation of the growth curves of Soderblom et al., 
1993). 

In Table lAT3l we report the presence of Ha and Nai(Di,D2) 
emission, and of a measurable lithium line at 6707.8 A. 
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Table A.l. Observing Log for the NEP spectra. Scan numbers are from Gioia et al. 2003. 



Scan 


Date 


Ra 


Dec 


Spectral type 


U (s) 


S/N 






(h m s) 










*20 


2004 Jun 28 


17 52 44.8 


+67 00 20 


G8 


5400 


70 


260 


2004 Jun 01 


17 58 28.3 


+67 26 08 


G2 


2400 


50 


330 


2004 Jun 02 


17 59 23.6 


+66 02 55 


G9 


1800 


70 


370 


2004 May 25 


18 00 02.2 


+66 45 53 


F8 


2400 


90 


*1420 


2004 Mar 08 


17 19 28.8 


+65 22 29 


KO 


3600 


60 


*1441 


2004 Jan 31 


17 20 05.6 


+62 06 22 


G2 


1200 


100 


1470 


2004 Mar 09 


17 20 26.7 


+67 03 37 


G8 


1800 


80 


1500 


2004 Apr 07 


17 21 10.9 


+69 48 02 


G2 


600 


200 


1511 


2004 Feb 03 


17 21 42.5 


+62 00 32 


K4 


3600 


60 


1580 


2004 Apr 07 


17 24 00.4 


+69 40 30 


K2 


3600 


50 


*1600 


2004 Mar 08 


17 24 26.8 


+64 12 23 


G2 


1800 


50 


1601 


2004 Mar 09 


17 24 38.8 


+64 40 51 


GO 


1800 


80 


1690 


2004 Mar 09 


17 26 45.4 


+69 37 53 


G2 


2400 


60 


1800 


2004 May 06 


17 29 39.4 


+68 47 38 


F5 


1200 


110 


1840 


2004 May 25 


17 30 19.9 


+69 55 27 


GO 


1800 


90 


1960 


2004 Apr 07 


17 33 16.9 


+67 12 08 


F2 


600 


170 


2050 


2004 Mar 09 


17 36 14.1 


+65 02 27 


KO 


2400 


70 


2100 


2004 May 06 


17 36 26.6 


+68 20 37 


M3 


600 


170 


2210 


2004 Jun 02 


17 39 16.1 


+70 20 09 


F8 


1200 


80 


2250 


2004 Mar 08 


17 39 56.1 


+65 00 04 


KO 


2400 


200 


2360 


2004 Jun 01 


17 42 26.5 


+69 07 58 


F8 


1200 


120 


2400 


2004 Mar 09 


17 43 01.6 


+66 06 46 


G2 


1800 


90 


2470 


2004 Jun 02 


17 43 51.7 


+70 31 39 


F3 


1800 


40 


2480 


2004 Jun 01 


17 44 00.6 


+70 15 27 


KO 


2400 


50 


*2970 


2004 Jan 31 


17 49 03.9 


+62 47 48 


F2 


600 


290 


*3080 


2004 Jun 28 


17 50 25.3 


+70 45 36 


K2.5 


1800 


90 


3220 


2004 Jun 01 


17 52 56.0 


+66 25 10 


KO 


600 


150 


3560 


2004 Jun 28 


17 57 03.7 


+68 49 14 


G8 


1800 


110 


3610 


2004 May 06 


17 58 01.4 


+64 09 34 


K4 


2700 


70 


3690 


2004 May 06 


17 58 48.0 


+63 50 39 


A2 


1200 


120 


3710 


2004 Jun 30 


17 59 13.8 


+64 08 33 


F2 


1200 


200 


*4090 


2004 Jun 28 


18 05 30.5 


+69 45 17 


G2 


2400 


50 


4100 


2004 May 25 


18 05 30.3 


+62 19 03 


KO 


1200 


150 


4130 


2004 Jun 02 


18 05 44.8 


+65 51 58 


F5 


1200 


80 


4180 


2004 May 25 


18 06 41.0 


+64 13 18 


KO 


600 


440 


4380 


2004 Jun 01 


18 08 45.4 


+62 56 37 


G5 


1200 


120 


4470 


2004 Jun 29 


18 09 55.8 


+69 40 39 


K3 


1200 


160 


4530 


2004 Jun 29 


18 10 49.9 


+70 16 09 


G9 


4800 


80 


4780 


2004 Jun 29 


18 13 48.6 


+68 31 32 


K2 


2400 


80 


4810 


2004 Jun 03 


18 13 51.4 


+64 23 57 


F5 


180 


300 


*4931 


2004 Jun 28 


18 16 21.2 


+65 29 39 


K8 


1200 


170 


4970 


2004 Jun 03 


18 16 49.7 


+65 04 26 


K2 


2400 


60 


5060 


2004 Jun 29 


18 18 31.9 


+70 42 17 


G5 


1800 


90 


5220 


2004 Jun 03 


18 20 19.2 


+65 19 19 


F8 


600 


120 


*5320 


2004 Jun 28 


18 21 46.8 


+63 57 10 


KO 


1800 


90 


5510 


2004 Jun 03 


18 25 09.0 


+64 50 21 


KO 


600 


260 


5520 


2004 Jun 03 


18 25 32.9 


+62 34 15 


K4 


2400 


60 


5950 


2004 Jun 30 


18 32 29.5 


+68 36 52 


GO 


600 


190 


6030 


2004 Jun 29 


18 33 47.8 


+65 13 33 


K3 


1200 


150 


6051 


2004 Jun 29 


18 34 33.6 


+69 31 45 


K5 


3600 


50 


6160 


2004 Jun 30 


18 36 13.5 


+65 29 15 


FO 


300 


230 


6163 


2004 Jun 29 


18 36 22.8 


+66 54 54 


G3 


600 


180 


6350 


2004 Jun 29 


18 39 25.4 


+69 02 54 


M3 


2400 


60 


6400 


2004 Jun 29 


18 40 44.4 


+70 38 47 


K2 


1200 


70 


6410 


2004 Jun 30 


1 8 40 56.6 


+62 44 54 


KO 


300 


300 


6510 


2004 Jun 29 


18 43 12.6 


+69 55 54 


IX) 


1200 


80 



*: We lack the Thorium- Argon lamp exposures for the nights in which these stars were acquired. They were wavelength calibrated using lamp 
exposures of the nearest night. 
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Table A.2. RESULTS: Vrad, Vrot, Li equivalent width, Li abundance, Binarity 
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nep 


Vrad (km s ') 


Vrot (km s"') 






Nu 




Spectral type 


notes 


*20 


-17.3±0.9 


<3.5 


<4.3 


5310+150 


<0.68 


<0.77 


G8 




260 


-19.5±1.1 


11.1+2.1 


112.4 


5790+150 


2.76+0.24 


2.71+0.19 


G2 




330 


-35.4±1.0 


6.7+0.8 


<3.6 


5250+150 


<0.53 




G9 




370 


-15.7±1.3 


7.0+1.0 


119.3 


6250+300 


3.22+0.28 


3.09+0.20 


F8 




*1420 




22.7+1.0 


- 


5150+150 






KO 




*1441 


-15.7±1.5 


4.0+1.4 


56.2 


5790+150 


2.34+0.20 


2.34+0.18 


G2 




1470 


-10.2±1.0 


7.4+0.4 


106.2 


5310+150 


2.22+0.24 


2.27+0.21 


G8 




1500 


-17.0±0.8 


5.5+0.7 


10.7 


5790+150 


1.57+0.18 


1.59+0.18 


G2 




1511 






87.0 


4450+350 


1.01+0.50 




K4(?) 


X 


1580 






65.1 


4830+300 


1.36+0.40 


1.55+0.37 


K2(?) 


X 


*1600 


-18.8±0.9 


5.4+1.4 


- 


5790+150 






G2 




1601 


-13.0±1.4 


11.0+0.7 


135.6 


5940+150 


3.05+0.24 


2.94+0.18 


GO 




1690 


-37.0±2.7 


10.1+1.0 


- 


5790+150 






G2 




1800 


-58.0±4.1 


16.3+1.4 


98.1 


6650+300 


3.28+0.22 


3.16+0.18 


F5 




1840 




>30.0 


- 


5940+150 






GO 




1960 




>50.0 


- 


7000+300 






F2 




2050 


-7.6±1.8 


12.8+1.0 


77.2 


5150+150 


1.85+0.23 


1.96+0.21 


KO 




2100 






- 


3500+350 






M3 




2210 


-7.1±1.8 


25.0+2.0 


- 


6250+300 






F8 (?) 


X 


2250 


-7.4±0.9 


3.8+0.7 


<3.0 


5150+150 


<0.33 




KO 




2360 


-26.2±1.4 


8.2+0.7 


44.9 


6250+300 


2.63+0.23 


2.58+0.16 


F8 




2400 


-14.3±0.9 


4.5+1.4 


<3.9 


5790+150 


<1.12 


<1.16 


G2 




2470 




8.3+1.0 


- 


6950+300 






F3 




2480 




>40.0 


- 


5150+150 






KO 




*2970 




>60.0 


- 


7000+300 






F2(?) 


X 


*3080 


-12.9±1.7 


20.6+1.0 


147.6 


4800+350 


1.85+0.49 


2.00+0.53 


K2.5 (?) 


X 


3220 


-25.1±1.0 


<3.5 


<5.4 


5150+150 


<0.60 


<0.72 


KO 




3560 


-2.6±1.3 


7.9+1.0 


124.9 


5310+150 


2.34+0.25 


2.38+0.24 


G8 




3610 


+9.7±1.0 


5.1+0.7 


13.9 


4450+400 


0.11+0.56 




K4 




3690 






- 


9000+800 






A2 




3710 




>30.0 


- 


7000+300 






F2 




*4090 


-40.7±2.2 


<3.5 


- 


5790+150 






G2 




4100 


+10.0±0.9 


4.5+0.7 


18.2 


5150+150 


1.14+0.18 


1.27+0.31 


KO 




4130 




>40.0 


- 


6650+300 






F5 




4180 


-13.6+1.0 


<3.5 


<3.0 


5150+150 


<0.33 




KO 




4380 


-21.8±1.1 


6.3+1.4 


29.0 


5560+200 


1.81+0.24 


1.86+0.23 


G5 




4470 


-13.5±0.8 


5.3+0.8 


24.3 


4800+350 


0.85+0.45 


1.06+0.44 


K3 




4530 


-21.6±1.8 


14.2+0.7 


228.7 


5300+150 


2.97+0.40 


2.84+0.26 


G9 




4780 


-26.3±2.1 


11.7+1.1 


422.3 


4830+300 


3.60+0.75 


3.40+0.65 


K2 




4810 


-37.0±1.6 


7.8+1.0 


- 


6650+300 






F5 




*4931 


+2.5±1.7 


3.6+1.0 


<5.5 


4000+150 


<-0.8 




K8 




4970 


+ 14.3±0.9 


<3.5 


<8.6 


4830+300 


<0.41 




K2 




5060 


-2.3+0.9 


7.6+2.1 


97.9 


5560+200 


2.44+0.26 


2.45+0.22 


G5 




5220 


-18.9±2.0 


25.8+2.2 


<4.8 


6250+300 


<1.61 


<1.59 


F8 




*5320 


-19.4+0.8 


4.0+0.7 


38.5 


5150+150 


1.49+0.21 


1.62+0.20 


KO 




5510 


-53.5+3.0 


20.4+2.1 


371.3 


5150+150 


3.79+0.67 


3.46+0.54 


K0(?) 


X 


5520" 


+8.0+1.3 


7.6+1.0 


93.2 


4450+350 


1.05+0.51 




K4 (?) 


SB2 


5520' 


+25.7+2.3 


4.2+1.4 


39.8 


4450+350 


0.60+0.49 




K4(?) 




5950 


-20.6+1.4 


11.4+1.0 


<4.3 


5940+150 


<1.30 


<1.31 


GO 




6030 


-13.2+1.0 


4.5+1.4 


<5.0 


4780+350 


<0.11 




K3 




6051 


-9.1+1.3 


15.1+1.0 




4410+400 






K5 




6160 








7300+300 






FO 




6163 


-7.4+0.5 


3.7+1.4 


49.4 


5700+150 


2.19+0.19 


2.21+0.18 


G3 




6350 




<3.5 




3500+350 






M3 




6400 


-23.9+0.8 


<3.5 


<5.0 


4830+300 


<0.17 




K2 




6410 


-25.6+1.2 


<3.5 


10.5 


5150+150 


0.89+0.21 


1.02+0.21 


KO 




6510 








730()+3()() 






I'O 





NOTE: An x in column eight indicates that the star is (or it is supposed to be) binary. We indicate SB2 for double-Uned spectroscopic binary. 
*: We lack the Thorium- Argon lamp exposures for the nights in which these stars were acquired. They were wavelength calibrated using lamp 

p.xnnsiire.s nf the. ne.are.sf nipht. 
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Table A.3. Presence of Ha and Nai(Z)i,D2) emission and presence of lithium feature (6707.8 A) 



nap binary Ha Nal(Di,D2) Lil 

*20 - - ~ 

260 - - X 

330 ... 

370 - - X 

*1420 X X 

*1441 - - X 

1470 - - X 

1500 - X X 

1511 X - X X 

1580 XX X X 

*1600 X - X 

1601 - - X 

1690 - X 

1800 - X X 
1840 
1960 

2050 - X X 
2100 

2210 X - 
2250 

2360 - - X 
2400 

2470 - X 

2480 X X 
*2970 

*3080 XX - X 

3220 

3560 - X X 

3610 - X X 
3690 

3710 - X 

*4090 X X 

4100 - - X 
4130 
4180 

4380 - - X 

4470 - - X 

4530 - - X 

4780 - X X 

4810 
*4931 

4970 - X 

5060 - - X 

5220 

*5320 - - X 

5510 XX X X 

5520 SB2 - X X 
5950 
6030 

6051 X X 
6160 

6163 - - X 

6350 X X 

6400 

6410 - - X 
6510 



*: We lack the Thorium-Argon lamp exposures for the nights in which these stars were acquired. They were wavelength calibrated using lamp 
exposures of the nearest night. 



